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Modeling scalar flux and the 
energy and dissipation equations 

By A. Yoshizawa 1 


Closure models derived from the Two-Scale- Direct-Interaction Approximation have 
been compared with data from direct simulations of turbulence. In the working 
session, we have restricted our attention to 1) anisotropic scalar diffusion models, 
2) models for the energy dissipation equation, and 3) models for energy diffusion. 


1. Anisotropic eddy-diffusivity model for turblent scalar flux 

The scalar flux is represented by a gradient diffusion model 


u\e> = -D i} 
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with a diffusivity tensor Dij that depends on the mean strain and vorticity tensors 
(Yoshizawa, 1985). 
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The accuracy of the model is shown in table 1 by comparison at several times t 
of the actual fluxes with the modeled fluxes in a direct numerical simulation of 
homogeneous turbulence in uniform shear S having a uniform scalar gradient. The 
scalar diffusivities D 22 , Dn, and D 21 are represented well by the model but D 33 
and flu are not. The performance of the model might be improved by the inclusion 
of unsteady terms in Ck suggested by the TSDIA analysis. 
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Table 1. Evaluation of the Scalar Diffusion Model from case C128U 
of Rogers, Moin, and Reynolds (1986) 
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.068 

.090 

1.97 

2.38 

1.32 

1.20 

.737 

1.98 

.936 

5.58 

10 

.102 

.108 

2.74 
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1.27 

1.23 
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1.98 
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6.52 

12 

.150 
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1.24 
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1.82 

.760 

6.60 

14 

.205 
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1.77 
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6.44 

16 

.250 
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1.66 
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2. A model for the dissipation of kinetic energy 

Here we contrast the familiar k — c model (model 1), 

£-c.,ir-c£ + A (C „ 

C (1 ~ 1.4, C' t2 ~ 1.9, C ({ ~ 0.07, 

with the model derived via the TSDIA approach (model 2; Yoshizawa, 1987) , 

w,= c 4 p - c 4+ c -' k f* t £ )2+ — ■ 

C ( 1 = C',2 ~ 1.70. 

Note that for both models, the eddy-viscosity approximation implies that 

f du4 dui s o 

c„-p = c, iC .vm~ + ^ , 

where C7 V ~ 0.09. For homogeneous turbulence in uniform shear 5, model 1 reduces 
to 


de - r ‘p r ‘ 

Wt~ c "k P ~ c ' 2 T ’ 


and model 2 reduces to 


% =c <'i p - c *T + 2C '" ks2 ■ 

These two models were tested against the homogenous shear turbulence fields of 
Rogers et al (1986) for 8 < St < 14. The resulting “constants” were found to be 
.97 < C ( i < 1.2 for model 1, and 1.7 < C’ f , < 1.9, -.025 < < -.018 for model 2. 

Note that the negative value of C' fi implies that the effect of rotation, given by the 
third term of the model, acts to reduce the dissipation rate. The simulation data 
also support the relationship C ( i = C \ 2 — 1.7 suggested by TSDIA. 

3. A model for the diffusion of kinetic energy 

The diffusion term 

Dk = 

in the equation for kinetic energy is usually modeled as (model 1) 

a v dk 
Dk ~ -d^ [CK Td^> 
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FIGURE 1. The eddy viscosity distribution in turbulent channel flow. Ex- 
perimental data of Hussain & Reynolds: A R = 32300, o R = 13800. Model: 
R = 32454, R = 12581, R - 3666. 


whereas the TSDIA analysis (Yoshizawa, 1982) indicates the presence of a cross- 
diffusion term 
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These models have been compared with the turbulent channel flow data of Kim et al 
(1987) (hereafter KMM) for 100 < < 180. The Reynolds number is 3300, based 

on channel half-height and centerline velocity, and the centerline is at y + — 180. 

When the data was fit with a single term of the model, the constants were esti- 
mated to be .11 < Ckk < -12 when C'ki was set to zero, and .06 < Ck* < -08 when 
Ckk was set to zero. At high Reynolds number, the eddy viscosity distribution 
(ut = u'v'/S) has maxima off the centerline of the channel. The cross-gradient 
term in model 2, when incorporated into a k — e model, can produce the off-axis 
maxima whereas model 1 cannot. However, at the low Reynolds number of the 
simulation, the eddy viscosity did not exhibit the off-centerline maxima strongly 
enough to allow the two constants to be found simultaneously from the data alone. 

If the constants are taken as Ckk = 08 and Ck ( = .03, the locations of the 
maxima and their values are reproduced. The data of KMM indicate a maximum 
ut/u = 16 at y/d = ±.5 while the model gives a maximum of 18 at y/d = ±.47. A 
comparison of the eddy- viscosity distribution of model 2 and experimental data at 
higher Reynolds numbers is shown in figure 1. 
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